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Abstract. — A  description  of  a  FORTRAN  IV  computer  program  for 
performing  sequential  analysis  on  four  common  distributions  after  the 
underlying  probability  distribution  is  known. 


Decision-making  and  cost-reduction  are  key  concepts  in  management 
research.  Sequential  hypothesis-testing  procedures  that  were  developed 
for  quality  control  in  industrial  processes  (2,  3,  5)  have  been  applied 
in  a  limited  way  to  biological  problems.  Forest  entomologists  have  used 
the  sequential  probability  ratio  test  of  Sobel  and  Wald  (2)  to  classify 
insect  populations  (4)  and  to  evaluate  the  need  for  insect  control  (l). 
This  note  describes  briefly  a  computer  program  that  was  developed  to 
eliminate  much  of  the  computational  work  associated  with  sequential 
analysis. 

Sequential  analysis  is  popular  with  forest  entomologists  because,  on 
the  average,  it  requires  fewer  observations  for  fixed  type-1  («)  and 
type-2  (P)  errors  than  fixed-sample  size  plans  require.  The  smaller  the 
difference  to  be  detected  at  a  given  type-1  («)  and  type-2  (P)  error,  the 
larger  the  average  sample  size.  In  this  method  of  analysis  we  assume 
that  the  observations  are  independent  and  the  probability  distribution  is 
known  beforehand. 

To  perform  any  sequential  analysis,  the  following  must  be  known: 

1.  A  description  of  the  underlying  probability  distribution. 

2.  A  statement  of  the  hypothesis  that  is  to  be  tested. 

3.  The  type-1  (a)  and  type-2  (/?)  errors. 

The  computer  program  was  developed  for  the  three-decision  problem 
(light,  medium,  heavy)  that  is  commonly  encountered  by  forest  ento- 
mologists. By  repeating  selected  input  data  in  certain  columns  of  the 


second  control  card,  the  program  can  be  used  to  obtain  the  required 
statistics  for  the  two-decision  problem  (control  vs  no  control). 

The  FORTRAN  IV  program  was  developed  by  the  authors  on  a  7094 
computer.  The  program  will  provide  the  following  statistics  for  the 
binomial,  negative  binomial,  poisson,  and  normal  distributions: 

1.  Equations  for  the  decision  lines. 

2.  Table  of  decision  boundaries. 

3.  Graph  of  decision  lines. 

4.  List  of  selected  points  for  the  operating-characteristics  (OC)  curve. 

5.  List  of  selected  points  for  the  average-sample-number  (ASN)  curve. 

For  the  normal  and  negative  binomial  distributions,  it  is  assumed  that 
the  variance  and  k-value  are  known  beforehand. 

Description  of  Control  Deck 

The  user  must  supply  the  following  control  cards  when  using  the 
program.  Two  control  cards  are  required  for  each  analysis  desired.  The 
first  control  card  contains  the  type-1  (a)  and  type- 2  errors,  the 

distribution  index,  and  the  number  of  decision  points  to  be  printed. 

Card  1        Column  Format 

1  -    3  F3.2       Type-1  error,  A 

4  -    6  F3.2       Type-2  error,  B 

7-8  12  Distribution  index 

01  -  Binomial 

02  -  Negative  binomial 

03  -  Poisson 

04  -  Normal 

9-11  13         Number  of  decision  points 

to  be  printed 

The  second  control  card  contains  the  parameters  needed  for  the 
analysis.  The  format  of  this  control  card  depends  upon  the  distribution. 
Binomial  distribution: 

Ho:  P<Pi 
Hi:  P2<P<P3 
H2:  P>P3 


Column  Format  Description 

1  -    8  F8.4  V1 

9-16  F8.4  P2 

17  -  24  F8.4  P3 

25  -  32  F8.4  P4 
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Negative  binomial: 

Ho:  M<Mi 
Hi:  M2<M<M3 
H2:  M>M4 


Column 
1  -  8 
9  -  16 
17  -  24 
25  -  32 
33  -  40 


Poisson: 

Ho:  A<X! 
Hi:  A25CA5C  A3 
H2:  A>A4 


Format 
F8.4 
F8.4 
F8.4 
F8.4 
F8.4 


Description 


Mr 

Ms 
K- 


" Index  of 
aggregation' 


Column 
1  -  8 
9  -  16 
17  -  24 
25  -  32 


Format 
F8.4 
F8.4 
F8.4 
F8.4 


Description 
Ai 
A2 
A3 
A4 


Normal: 

Ho:  p<h 
Hi:  ^2</a</x3 
H2:  /a>/a4 


Column 
1  -  8 
9  -  16 
17  -  24 
25  -  32 
33  -  42 


Format 
F8.4 
F8.4 
F8.4 
F8.4 
F10.4 


Description 
Pi 

1*2 
^3 

m 

a  -  standard  deviation 


Example 

Suppose  we  wish  to  calculate  the  sequential  plan  described  by  Waters 
(4)  for  the  spruce  budworm.  The  basic  data  for  this  plan  were  fitted  to 
the  negative  binomial  distribution.  Waters  listed  the  following  parame- 
ters for  his  plan: 


a  =  0.10 
Mi  =  3.00 
M3  =  9.00 


p  =  0.10 
M2  =  6.00 
M4  =  12.00 


7.288 
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Figure  1.— Program  output  of  decision  equations  and  se- 
quential table. 
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Figure  2.— Graph  of  decision  line  equations. 
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TABLE  OF  VALUES 

FOR  THE  OC 

-AN_°_ ASH  CURVES—  NEG  BIN 



CASE  1 
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CASE  2  
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LCk       VS  MEDIUM 

MEOIUM 
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OC 
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P 

OC 
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0.486 

3.730 

1.432 
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MEAN   1  -     3.0000  CC-VALUE  «  0.900  MEAN  2  -     6.0000  OC-VALUF  *  0.100 

MEAN  3  »     9.0000  CC-VALUE  «  0.900  MEAN  4  «   12.0000  OC-VALUE  »  0.100 

Figure  3.— Program  output  of  OC  and  ASN  values. 


From  this  information  the  control  deck  would  be  prepared  as  follows: 

COLUMN  NO. 
Card 

No.   I  2  3  4  5  6  7  8  91  0  1  2  3  4  5  6  7  8  92  0  1  2  3  4  5  6  7  8  93  0  1  2  3  4  5  6  7  8  94  0 

1  01001002025 

2  000300000  0060000000  9000000120  0000007228  0 

Figures  1,  2,  and  3  show  the  sequential  plan  produced  by  this  program. 
The  decision-line  equations  produced  by  this  program  differ  slightly 
from  those  published  by  Waters,  but  these  differences  are  probably  due 
to  rounding  errors.  The  occasional  discrepancies  between  values  in  the 
sequential  tables  are  also  the  results  of  rounding  errors. 

A  program  deck  with  examples  can  be  obtained  from  the  authors. 
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